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ABSTRACT 



A pertubation expansion is formulated for the three- 
dimensional, nonlinear, acoustic-wave equation with dissipa- 
tive term describing the viscous and thermal energy losses 
encountered in a cylindrical cavity. The theoretical results 
show that nonlinear effects in sloshing modes are strongly 
suppressed . 
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c 

Co 

c„ 

Cf 

-k 



M 



p. 



u 

C>C 

C^TlCo 

f-e 

Uj 

l/U^ 



phase velocity in the cavity 

-Cxf - a'/at 

‘<i e-e- 

effective phase speed associated with a standing wave 
at resonance 

specific heat at constant pressure 

propagation constant associated with a standing 

wave 

peak Mach number of the driven standing wave 
= P - p) = - “ acoustic pressure 

instantaneous and equilibrium total pressure in the 
field 

absolute temperature in kelvin 
particle velocity 

particle speed of the 7]'th_order perturbation solution 
infinitesimal-amplitude attenuation constant 
temporal decay constant of a resonance 
=(T + 1) /2 for a gas 

=^p/Cv 5 ratio of specific heats for a gas 
instantaneous and equilibrium densities 
velocity potential 

=zzf, (angular) grequency at which the cavity is driven 
(angular) frequency of a resonance 
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I. INTRODUCTION 



The topic of finite amplitude acoustic standing waves in 
sloshing modes of a cylindrical cavity is interesting 
theoretically, but development of the subject has not been 
extensive. The purpose of this research is to present the 
results of a power series perturbation approach to the 
problem . 
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II . THE NON-LINEAR WAVE EQUATION 



A. GENERAL 

It is well known [l] that for <^k^and M«l, where 
measures the fractional loss per wavelength and M is the peak 
Mach number of the source, loss terms and nonlinear terms in 
the constitutive equations can be separately approximated with 
the help of linear, lossless acoustic relations. The force 
equation appropriate for acoustical processes in systems for 
which gravitational effects are unimportant is 

+(U'V)u + -^yp (2-1) 

where (_J = LI is the particle velocity, ^ is the instantaneous 
density of the fluid, P is the instantaneous total pressure in 
the field, and the operator symbolically describes those 
physical processes leading to absorption and dispersion. We 
used two additional equations.' The first is the equation of 
state for a perfect gas 

P - ^TTk (2-1-a) 

where r is a constant whose value depends on the particular 
gas involved and is the absolute temperature in Kelvin. 

The second is the continuity equation 
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( 2-1-b) 



IL 



+ V- [U + s)U ] = 



where s = ( is the condensation at any point and ^ 

is the equilibrium density of the fluid. If we ignore 
rotational effects, then 



a =Vi 



( 2 - 2 ) 



where ^ is the velocity potential. Combination of Eqs.(2-D- 
(2-2) and the neglect of terms of orders higher then , M(°^/-|;) , 
and ( yields a quadrat ically nonlinear wave equation, 



D I o' r y ^ U ^ 

C It = -1 I + ^12^ 









(2-3) 



where Co^ = (^|/^p) (adiabatic), p=acoustic pressure , 

is the ratio of heat capacities, and 







(2-4) 



The left-hand side of Eq.(2-3) is the classical, linear wave 
equation with losses pertinent to the system under study. 
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The right-hand side can be interpreted as a forcing function 
consisting of a three-dimensional spatial distribution of phase- 
coherent sources. 

In a second-order perturbation theory, this volume forcing 
function is obtained from the classicaK first-order) solution 
P]_ of the acoustical problem. The second-order perturbation 
solution P 1 +P 2 describes the non-linearities P 2 resulting from 
the self interaction of the classical solution P]_. 

Higher-order perturbation solutions consider the interaction 
of the nonlinear solution with itself, and the forcing function 
is composed of products of both classical and non-linearly 
generated terms. Thus, if a system is driven at frequency w, 
the non-linear term in Eq.(2-3) will force the existence of 
all integer multiples Tllli' of the driving frequeney and the full 
solution must contain all harmonics of the input frequency. 

In a closed cavity, each of those nonlinearly generated waves 
whose frequency lies near the resonance frequency of a standing 
wave of the cavity and whose assciated spatial function matches 
that of the standing wave can be strongly excited [ 2 ] . 

As far as the author has been able to determine, there has 
been only one previous study of this system published in the 
open literature. This was by Maslen and Moore in 1956 [s] . 
Their approach resulted in a series expansion which did not 
converge if the relevant nomal mode frequencies were integerally 
related. Their interpretation of the quenching of the nonlinear 
effect was based on the "scattering effect of the wall". 
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We feel our interpretation based upon nonlinearally generated 
volume sources stimulating the allowed standing waves of the 
cavity is more accessible and informative. Further, the 
mathematical approach developed herein appear to avoid any 
difficulties in convergence. Their conclusion that high 
amplitude monofrequency transverse oscillations can exist is 
consistent with our finding. 

B. APPLICATION TO THE CYLINDRICAL CAVITY 

The circular cylinder is provided with a "point" source of 
sound. By properly positioning the cylinder with respect to 
the sound source it is possible to effectively drive the 
enclosed air into various modes of vibration. The rigid- 
walled cylinder forces the component of particle velocity 
perpendicular to each cavity surface to vanish at the surface. 
The resulting steady-state solution to the linear wave equation 
in cylindrical coordinates is 

(2-5) 



where Jj-^ are the cylindrical Bessel functions and application 
of the boundary condition to the sides yields 




Jn-m 

a. 



(2-5-a) 
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where a is the radius of cylinder and are the arguments 

.-til 



of the exterema of the 'Jj Bessel function. 

The nomal mode frequency is dependent on a.nd ^ 



rrnri 






firm 



-i -1C 






- 



(2-5-b) 



The standing wave well he identified by the ordered integers 
(n, m, X ) describing its spatial dependence. 

1 . Symmetric Modes 

The simplest waves in cylindrical coordinates are those 
that depend only on the distance r from Z-axis, the gradient 
takes the form 

where 

1 =^^ as-.: 



and k=w/c is the wave number or propagation constant. 
The Laplacian becomes 






(2-6-b) 



The suitable trial solution appropriate for sy.mmetric modes is 
the linear solutionCO, m, 0) of frequency w 
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The velocity potential for this limit can be defined as 




( ^ ) M Jo Cjf) 



( 2 - 8 ) 



and 



a 



/C is 



Ul = ^ 

Co 0 



ri ) CCS lut = - f M J, CC6 wt 



(2-9) 



Because we assume that the surfaces of the cavity are rigid, 
then at r=a the appropriate boundary condition is Jg'=0. 

Thus , we have -la = V 

otn • 

To generate the second order solution, we first note that 






( 2 - 10 ) 



and 




rrj" ( I -f COSZUJ^t) 



( 2 - 11 ) 



Next, the second derivatives of Eqs.(2-10) and (2-11) with 
respect to time are 



(2-13) 



2l r — 

^ Cc 






Substituting Eqs.(2-2) and (2-13) into the first term of the 
right-hand side of Eq.(2-3) yields 



(;cw-) 

4 



Nyi 



^ ^ —Z^ 

CoS Xtl' t , J, - Y Jo 



( 2 - 



and the second term of right-hand side of Eq.(2-3) becomes 

x'+ 4-ji'- 3r ) -( 



4 



(2. 



since 



c:v‘j/= cr-JoX 



( 2 - 



or 



X . ^ I “P ^ 1 I "-T“X Z 

Co V Ji Ct X T xTo - J, ) 



( 2 - 



14) 



15) 



16) 



17) 



With the use of Eqs.(2-14) and (2-15) we can write the right- 



hand side of Eq.(2-3) as 




•COS Zl^t 



(2-18) 



and the appropriate inhomogeneous wave equation for the second 
order perturbation solution is 







(2-19) 



2 . Non- Symmetric Modes 

The non planar waves in cylindrical coordinates are 
those that depend on the distance r from Z axis and the angle 
from X-axis. The gradient takes the form 




14 



and the Laplacian becomes 




- 4,^ C A- 



-h 



1 



? 



H- 






0 



(2-20-a) 



The suitable trial solution appropriate to the forcing function 
of frequency w for exciting the (n, m, 0) standing wave is 



R 







m 



lOt 




( 2 - 21 ) 



where is the Bessel function of the first kind and 

order n. The velocity potential is approximated by 




= ( -^ ) M 31 CCS wt 



( 2 - 22 ) 



and 



U./^^ is 

^ M C ^ cc« U't - q ^ X coiu>t) 



( 2 - 23 ) 



Thus , 







( 2 - 24 ) 



or 




( 2 - 25 ) 
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and 



+ [vf ( i + CC6 m ( J-Tii-i + J T)-i ) J^nX-i 

Co T 



(2-26) 



and the second derivative of Eqs.(2-25) and (2-26) with respect 
to time is 




(?W)_r 

i 



v]"j^"(l+dc6a'n5>)coS 




(2-27) 



and 






_ C-^'f jv|^ 

■ i ' 



( 



'on 



+ ( ^ 



-T) 






+ 





(2-28) 



Substituting Eqs.(2-26) and (2-27) into the first term of 



Eq . ( 2-3 ) yields 



(2curf 

g 




cos 





(2-29) 
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f 



or 



s 



;iwt 1 C Xn+I + J-IH ) 






'Vi 






V 




s ince 



and 






In order to solve the second term of Eq.C2-3), we can 

c/V^Jn" Jtr X > 

CoV W -tm «6 ^"’5 = xw"(X,3,h -x: 

-T^^X-^+ )C0S x-n^ 



(2-30) 



(2-31) 



(2-32 ) 



use 

(2-33) 



(2-34) 



(2-35) 
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since 





C37ltl 



(2-36) 



and 

(^) Jm “ " Jnti JdTj 



(2-37) 



Thus, Substituting Eqs.(2-33) through( 2- 35 ) into the second 
term of Eq.(2-3) yields 









.-t ) (4 J, 



Titl 



4- 

I 





(2-38) 



and 



— g" v"C 1 ( 1 + -«>t ) 

Z Oo -t 



J-/J-T 

X V ^ Jny 



lOHX 






- 3^ “ Xi-t-) i 3^41 Tt}-! ) ^ '3'ti 



f 






(2-39) 
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With Eqs.C2-30), (2-38) and (2-39) we can get the inhomogeneous 
wave equation for the second order perturbation as 














COS 




I -r- . I ^ ^ 

15 ^ Jq)^) 



If we let n=0 in Eq(2-40), it reduces to Eq.(2-19), as it must. 
Substituting n=l into Eq.(2-40) yields the equation with 
forcing term resulting from the (1, m, 0) sloshing mode, 







.2 




7 ^ 










ctS ;t il>-t • 




r 



^ 3a. Ct:: j^)3^"t'^3o 



; 





-L_\T^ 

g 




y 





(2-41) 
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III. METHOD OF SOLUTION 



Recall that the equation with forcing term resulting from 
thed, m, 0) sloshing mode can be written as a function of 
frequency , 



4 



■ jr+ y (4 “ II )J 



r 



l-rT-lTT-fX-t-'lT 

S J 5 J 1 ^JxJo U g 



CCS ^ 



4 



( 3-1) 



and the left-hand side of Ea.(2-3) is 



r - X 0 i ± L J — ^ \ 
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^ r 









R 



,Co' 



( 3-1-a) 



If the harmonics of the frequency at which the cavity is driven 
are not close to any resonant frequencies, we can ignore the 
lossy term. So, Eq. (3-1-a) can be written as 



Co 4 ^ 3 -^ ) 






1. 



(3-2) 



Let us assume that the solution of Eq.(3-1) can be expressed as 

P _ 






- = f3|^CCS ^4 



( 3-3) 
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and define 





M^Vn COS Xli’t 



and 



£iil = ri^ u-« c-^v) 



(3-4) 



( 3-5) 



Combination of Eqs.C3-2) and (3-3) yields 











~i — ^ T 

‘ f J4)H 




J 






J 



(3-6) 



and 





T T 

(g JiHciii-L 




( 3-7) 
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A. POWER SERIES METHOD 

In order to solve the right-hand side of Eq.(3-7), we 
use the definition [^ ] 

(^1/ ""z- 

4 i 

(-1 X r j 

p (V+'fe't'i ) n ) 1 ~k ! 

where P( I j = J j 
1 . Coefficient 

a. Substituting V=4i= n+2 into Eq.(3-8) yields 
— ^ M Z7)t4-^ 

■"-fe 

/ ! C^-r'fe + 4)/ "k? 

or 

X.(1)=z 

•4^ 



where 



n /■ 1 


( f 4- ) / 




^ N 

^ > 


( 'fe/ 



can 



(3-8) 



(3-9) 



( 3-10) 



( 3-11) 
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b. Substituting V= n+2 , /J_~n-2 into Eq.(3-8) yields 






/>o 

5 Z 

,2^ 






(3-12) 



or 



Co 



4:=o 






(3-13) 



where 



U = (-0 - 



-k c ^ ) - 






( 3-14) 



c. Substituting i) = n+l , .{( =n-l into Eq.(3-8) yields 



, I x^")' 

j 1HI ^ X-i ~ 1 ) 



CO 




-1=0 



n 

! (XX 



(dOti-l)/ 



(3-15) 



or 



Xti (1 )j-,i (1 ) = H (^1 ) 

-H 






(3-16) 
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where 






( 2.Q0 + J-fe ) 






or 



where 



d. Substituting ])= ^ =n into Eq.(3-8) yields 



z-k 



XcVHiifEj 



':?) I 



1-0 



_ { I \-fe -t-Ji-'fe '> / 



C‘ ' ) 






e. Substituting V =^(. =n+l into Eq.(3-8) yields 







X+, (X' = (,^ 

(^-l)^(x-Tid-z-fe-rx) / 

(nn-fe-M)! (xm-r-^xx) ! -fe ! 



( 3-17) 



(3-18) 



(3-19) 



(3-20) 



( 3-21) 
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or 



■j;, 

■^--0 






(3-22) 



where 






^ (j.Ol-t-2.-fe.T-S^ ) ( 




I kl 



(3-23) 



f. Substituting V ='</t=n-l into Eq.(3-8) yields 



where 



Xh(1) = (t1) Z- 



<>o 



(3-24) 



or 



CO 



H/rvfe (iV 

-h) 



^Cv-f4:-i ) 



( 3-25) 






_ / ,vfe ( -;c) ! 




( 3-26) 
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g. Substituting ]) -/it -n-2 into Eq.(3-8) yields 






or 



■^=0 




where 



0)i« ~ ) 



i 



( 2.'T)-t 2 k - i 






2. Vn(y) 



With the help of Eqs.(3-9) throughC 3-2 9 ) Eq.(3-7) 
be written as 




( 3-27) 



(3-28) 



(3-29) 



can 



(3-30) 
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I 




where 



A = 



X 



( 3-30-a) 



and 



a = 



CziL 



( 3-30-b) 



Now, let us assume that 



0 ° 



■V 



^ 1 ) 4 ( ^ 






( 3-31) 



The first derivative with respect to y= -fer can be expressed as 

Co 



dV>i I 4^ A /I 

j . ~ ^ ^ ( 01 +^) 

^ ^=0 



(3-32) 



and the second derivative with respect to y="^r can be expressed 
as 



iC'Df'kyz 



( 3-33) 



Substituting Eqs.(3-31) through ( 3- 3 3 ) into the left-hand side 
of Eq.(3-30) yields 
CO 




'k=o 



( 'kf j('D+t)('>1+42- y) Clfe j 
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-4 Q« Tt ) 






N 



y 




(3-34) 



Eq.(3-34) must be equal the right-hand side of Eq.(3-30). Thus, 

U,p-^AE J 

(3-35) 

where k=0,l,2,3 (3-36) 

If Eq.(3-35) is substituted into Eq.(3-36), then Eq.(3-35) can 
be expressed as 




4-.10 + I ) < I ) 

(u t(^f( ^'^+x)c 

c 

m 

t (.i^J c^'i+-fe+i )(-feti ) 




(X,,.hr ^ 

(3-37) 
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The normal component of particle velocity is equal to zero on 
the boundaries. From the application of the boundary condition 
at the rigid boundary at r=a, 

60 



d 



d ^ 



I 

4 



(zj ) =0 



r-r 






where 



Substituting Eq.(3-36) into Eq.(3-38) yields 



(}y)+i)Q,)l (4j4 + 3 



(3-38) 






- O 



( 3-39) 



3. Un(y) 



With Eqs.(3-9) throughC 3-29 ) , Eq.(3-6) can be expressed as 

[(■*1X4 + 1^1+0 U-‘1> - 






o(qi-h4) 










4(9)44-i) 



(3-40) 
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I 



! 

c 






where, as before, A = '^/2-l/8 and a = 



v\; 

cS 



Now let us assume that 



/>o 






2CD-i-^-Z) 



( 3-41) 



The first derivative with respect to y=-|ir can be expressed as 

-TT = -T / (nn-'fe-x; hvk Cx3 j 

^1 



( 3-42) 



and the second derivative with respect to y=-^j^can be expressed 



as 






iU:^ = -i y ( 3 43) 

d f ' k 



Substituting Eqs.(3-41) throughC 3-4 3 ) into the left-hand side of 
Eq.(3-40) yields 



Txh 



^=0 



II • 1 

) bt-fe ^ 









( il) 






( 3-44) 



Eq.(3-44) must equal the right-hand side of Eq.(3-40). 
Thus , 
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k-fe+i + kfe = (--jj 5 ^,h 



++H 



71 .-^-1 16 



)^f\ 



(3-45) 



where \ - 1, 2, 3 

If Eq.(3-45) is substituted into Eq.(3-36), then Eq.(3-45) can 
be expressed as 



/'Jll “ GtIo 



+4 ko 

1?T» + (^)* C'^+l) W ]6 

k HmT i^zj U,r 



i ^ n 



(3-46) 
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1 






( 

I 

I 

II 



Application of the boundary condition at r=a yields 



cl 









=■4 

j' 



-=o 



where 



y(4<c)= clj/d L^<^) 

Substituting Eq.(3-36) into Eq.(3-47) yields 



b 1 / i 

l;/no ( "z j ) + ('Tl-i )p«|(-cj ; 






O-Dr 1 



( 3-47) 



4. - -|- ) 



= 0 



(3-48) 
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IV. CONCLUSIONS 



Recall that we can compute from Eqs.(3-37) through ( 3- 3 9 ) 
and Eqs.(3-14), (3-17), (3-20). Thus, with the use of Eq.(3-31), 
we can get . From Eqs.(3-46) through( 3-48 ) and Eqs.(3-ll), 



( 3-20), ( 3-23), (3-26), ( 3-29 ), we can compute . Thus, with 

the use of Eq.(3-14), we can get . Therefore we can compute 
Vn(y) and U^(y). Let us calculate the finite amplitude effects 
resulting from the nonlinear distortion of a forced radial mode. 
If we excite a(0, m, 0) mode and obtain the pressure at the 
circumference, then ^ - = J ^ . Given this value of , use 



of Eqs.(3-14), (3-17), (3-20) give the quantities C, D, E. 

From Eqs.(3-35) and (3-36) we can compute Q'r':^ . Thus, we can get 
V^Cy) from Eq.(3-31). Let us calculate the finite amplitude 
effects resulting from the nonlinear distortion of a forced 
sloshing mode. If we excite a(l, m, 0) mode and obtain the pre- 



quantities B. E, G, H, 0. From Eqs.(3-45) we can compute . 
Thus, we can get Uj.^(y) from Eq.(3-41). Substituting n = 0 at 
radial modes into Vj^(y) and U^^(y) yields Vg(y)=Ug(y). Therefore, 
with the use of Eq.(3-3), we can get the solution of Eq.(3-2). 

A. RADIAL MODES (0, m, 0) 



ssure at the circumference. 




Given this value 



of , use of Eqs.(3-ll), ( 3-20 ), ( 3-23), ( 3-26) 



(3-29) give the 
compute . 
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then 






so that 







B. SLOSHING MODESCl, m, 0) 




then 



N 




SO that 
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